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Abstract 

Full expressions of the B® — > tt + tt~ and B® — > n ir amplitudes, given in terms 
of matrix elements of operators of the effective weak Hamiltonian, are used to study 
the dependence of the relevant branching ratios on the different contributions. The 
uncertainty in the extraction of the weak phase a from the measurement of the time- 
dependent asymmetry in B® — > tt + -k~ decays is also analyzed. We find that, among 
several effects which may enhance the B^ — * vr 7r° branching ratio, the most important 
is due to "charming penguin" diagrams that have never been studied before. These 
diagrams easily increase BR(B® — > 7r°7r°) up to a value of 1-3 x 10 -6 . The same 
effect produces, however, a large error in the extraction of a from the measurement of 
the B® — ► it + tt~ time-dependent asymmetry. We show that it is possible to determine 
charming-penguin amplitudes from the experimental measurement of many decay rates. 
Their effect is impressive in B + — ► ir + K° and B® — ► K + ir~ decays, where charming- 
penguin contributions easily give values of BR{B + — > tt + K®) and BR(B® — > K + ir~) 
of about 1 x 10~ 5 . Among other possibilities, we also suggest to use B^ — > K°K°, the 
BR of which can be as large as 2-3 x 10 -6 , to determine the size of charming-penguin 
amplitudes. 



* On leave of absence from INFN, Sezione Sanita, V.le Regina Elena 299, Rome, Italy. 



1 Introduction 



The study of B® — > 7T7t (5^ — > 7T7r) decays is of paramount importance for our un- 
derstanding of CP-violation in the Standard Model and beyond. In particular the 
measurement of the time-dependent asymmetry 

jVQgg ^ vr+vr-)(t) - N(B° d -> 7r+7r-)(t) 

1 J JV(5° - 9r+7T-)(t) + iV(£° - TT+TT-Xi) 

(1 - |A| 2 ) cos(AM d t) - 2ImAsin(AM (i t) 

i + |AP ( > 

may allow the extraction of the CP-violating phase a, see for example [Q]. The cleanest 
method to extract sin 2a is from the measurement of the asymmetry, combined with 
the separate determination of the 1 = and 1 = 2 decay amplitudes, including the 
relative phase (2j. These can be obtained by measuring the B + — > 7r + 7r°, B® —* tt + tt~ 
and B® —* 7r 7r° (and the corresponding ones in the B® case) branching ratios. With 
these measurements, we get rid of our ignorance of the hadronic matrix elements of the 
weak Hamiltonian. Unfortunately, most of the theoretical analyses tend to predict a 
very small B^ — * branching ratio, thus making the model-independent extraction 
of sin 2a impossible in practice. 

If sin 2a has to be extracted from B^ — > 7r + 7r~ only, the main uncertainty comes 
from the contribution proportional to At = VtdV^, which is usually called "penguin 
pollution". In several studies, the decay rates and the uncertainty of sin 2a have been 
estimated by using some specific model to evaluate the hadronic matrix elements of 
the four-fermion operators entering the effective weak Hamiltonian ||— p|. In the most 
popular approaches the amplitudes have been computed by assuming the factorization 
hypothesis. The matrix elements of the weak currents necessary for the evaluation 
of the factorized amplitudes are then taken from a specific quark model or from the 
HQET §-[§. 

In this paper, we present a "model-independent" analysis of the uncertainty on 
sin 2a and of the ratio R = T(B% -> w°ir )/r(B% -> vr+vr-). By "model-independent" 
we mean that we do not make specific assumptions on the hadronic matrix elements 
of the operators, such as factorization or the absence of final state interactions (FSI). 
On the basis of simple "qualitative" physical considerations, we allow, instead, the 
matrix elements to vary within certain "reasonable" ranges, and check the stability of 
the results against such variations. This is particularly relevant for R, because of the 
delicate cancellations occurring between different amplitudes present in the B d —> 
case. Indeed, for this decay, the assumption of factorization and of the absence of FSI, 
or any approximation used to predict the value of the amplitude, may lead to an 
underestimate of the value of the decay rate. 

Our calculations are based on complete expressions of the decays amplitudes for 
B + — > 7r + 7r°, B® — > ir + ir~ and B® — > 7r°7r°, given in terms of diagrams representing 
Wick contractions of the operators of the effective Hamiltonian between the relevant 
external states. These formulae allow us to clarify assumptions and approximations 
usually made to evaluate the amplitudes, which have not been spelt explicitly in pre- 
vious studies. In particular, we show the presence of diagrams, involving operators 
containing charmed quarks (defined as Q\ and Qi in sec. that contribute to the 



penguin pollution, and that have never been considered before. We call these dia- 
grams "charming penguins" . 

Among several effects which are able to enhance the B® — ► 7r°7r° branching ratio, 
the most remarkable is due precisely to charming penguins. Their contribution may 
increase the estimate of BR(B® — > 7r°7r°) up to a value of 1-3 x 10~ 6 . The reason 
is that, unlike the case of the penguin operators Q3-Q10 which have small Wilson 
coefficients (of order a s /127rln(mj//i 2 )), the coefficients of Q\ and Q2 are of 0(1) and 
there is no reason to believe the corresponding matrix elements to be small. Charming 
penguins are also relevant for the B® — ► ir + ir~ amplitude and may give a large shift 
A = sin 2a e f ' —sin 2a ~ 0.4-0.8 between the physical value of sin 2a and the "effective" 
value, sin2a e -^, which can be extracted from the experimental measurement of ImX. 
As a comparison, when charming penguins are not included, the typical value is A ~ 
0.1. As A increases, however, also BR(B® — ► 7r°7r°) becomes larger, thus opening the 
possibility of extracting sin 2a with the isospin analysis proposed in ref. 0]. 

We finally show that many decay rates are expected to be dominated by charming- 
penguin diagrams. Among the various possibilities, we consider B® — ► K°K°, B + — » 
Tt + K° and B^ — > K + ir~ decays. In these cases, we give explicit formulae for the 
amplitudes, show that the largest contributions are those expected from charming 
penguins and estimate the corresponding branching ratios. 

The most impressive effect of charming penguins is found in B + — ► tt + K° and 
B® — > K + ir~ decays. Assuming reasonable values for the charming-penguin contri- 
butions, we find that their branching ratios may even become larger than BR(B® — ► 
7r + 7r~). This observation is particularly interesting because, in absence of charming- 
penguin diagrams, the B + — ► tt + K° and B® — ► K + it~ rates turn out to be very 
small either because there is a Cabibbo suppression or because the non-Cabibbo sup- 
pressed terms come from penguin operators which have rather small Wilson coeffi- 
cients (unless the corresponding matrix elements are exceedingly large). While fin- 
ishing this analysis, we were informed that the CLEO collaboration has measured 
BR{B° d -> K + k-) = (1.5i°;4 ± 0.2) x 10 ~ 5 II- The prediction that charming-penguin 
diagrams are important and give large B + — > ir + K° and B® — > K + ir~ decay rates is 
supported by this measurement. By using the experimental information, we predict 
BR(B + — > ir + K°) ~ 1 x 10 -5 and we call for a search of this decay mode. 

Other interesting decay channels, where charming penguins are expected to play 
an important role, such as B d — > ir°r] (or ?/), B® — > cpK° or B® — > r\r\ decays, will be 



extensively discussed elsewhere [10 



The plan of this paper is the following. In section we introduce the effective 
Hamiltonian given in terms of four-fermion operators, and of the corresponding Wilson 
coefficients; we also define the full set of diagrams in terms of which the B — * tttt ampli- 
tudes can be expressed; the final formulae of the different amplitudes are given at the 
end of this section. Formulae and approximations for the B® — ► K°K°, B + — > tt + K° 
and B® — ► K + tt~ amplitudes are discussed in section ||. In section ^ we present several 
physical arguments which are used to guide us in estimating the matrix elements; we 
also explain the main criteria used in the numerical analysis. In section ||, we give and 
discuss the main numerical results for R, sin 2a and for the B® — ► K°K°, B + — ► tt + K° 
and B^ — > K + ir~ branching ratios. 



2 Relevant formulae for B — > 7T7T decays 



The effective weak Hamiltonian relevant for B — > 7T7T decays is given by 

Htff 1 = X u ^[(c 1 ( f ,)(Q u l (^)-Qi(^)) + C 2 (^)(Q u 2 ( f ,)-Q 2 ( f ,)) 

+ rC^)-Q(p)} (2) 

where Q(fi) = (Qi(fi), Q 2 (v), ■ ■ • , Qio C?(aO = (Ci(/i), C 2 (/x), . . . , C 10 (//)), A u = 
V u dV* b and similarly A c and A^; r = —\t/\ u an d /i is the renormalization scale of 
the operators Qj. A convenient basis of operators [11|-[14], when QCD and QED 
corrections are taken into account, is 



Qi 


= (P d ){V-A){uu)(y_ A ) 


Q 2 


= (M)(y-A){ud){v-A) 


Qi 


= (bd)(y-A)(cc)(V-A) 


Q2 


= (bc)(y-A)(cd)(y-A) 


Q 3 ,5 


= (bd)(y-A)J2^)(v^A) 
q 




= J2(b<l)(V-A)(qd)(V-A) 
q 




= ~ 2 J2( b ~ < l)(s+P)(qd) (s _p ) 
q 


Q7,9 


3 

= ~( bd )(V-A) e q(v<l)(v±A) 
Z q 


Qs 


= -3j2e q (bq) {s+P) (qd) (s _p ) 
q 


Q10 


3 



(3) 



where the subscripts (V ±A) and (S±P) indicate the chiral structures and e q denotes 
the quark electric charge (e u = 2/3, = —1/3, etc.). The sum over the quarks q runs 
over the active flavours at the scale fi. 

Wick contractions of 7~C e ff between hadronic states give rise to the diagrams shown 
in figs. these are "Disconnected Emission" (DE), denoted also as T or T' in 

ref. |0|]; the colour suppressed (non-factorizable) "Connected Emission" (CE), de- 
noted also as C or C; "Disconnected Annihilation" (DA), denoted also as A or A'; 
"Connected Annihilation" (CA), denoted also as E or E'\ "Disconnected Penguin" 
(DP), "Connected Penguin" (CP), "Disconnected Penguin Annihilation" (DP A), 
"Connected Penguin Annihilation" (CPA). We assume SU(2) isospin symmetry. For 
penguin diagrams, we introduce a label which identifies the quark flavour in the pen- 
guin loop. Thus, for example, DP(s) denotes the disconnected penguin diagram DP 
of fig. |2], with a strange quark in the internal loop (a similar notation is adopted also 
for CP, DP A and CPA). Since we do not distinguish up and down quarks, we simply 
call DP (CP, DP A and CPA) those diagrams with an up or a down quark in the 
loop. In eq. (0), different Dirac structures appear, namely L = (V — A) x (V — A), 



DE CE 




DA CA 



Figure 1: Non-penguin diagrams. The dashed line represents the four-fermion operator. 




DP CP 




DPA CPA 



Figure 2: Penguin diagrams. 



R = (V - A) x (V + A) and S = (S + P) x (S - P) 1 . Thus, for example, the notation 
CEl and CE$ denote the connected emission diagrams with a (V — A) x (V — A) 
or (S + P) x (S — P) operator inserted, respectively. We are now ready to give the 
complete expressions for the decay amplitudes under study. It is convenient to write 
their AI = 3/2 and AI = 1/2 components separately 

A +0 = A(B+^^°) = ^\ U -^A 2 , (4) 

A+- = A(B d ^n + n-) = ^=\ u (A 2 -A ) , (5) 

A oo = i( B o^ A o ) = ^ Ai(( V2i 2 + -Li ) i ( 6 ) 



where 
and 



A 2 = A 2 l + rA 2 , A = A% + tAI, (7) 



A u 2 = ~[C 1 + C 2 ](DE L + CE L ) (8) 

A\ = -^[C 7 (DE R + CE R )-2C 8 (DE S + CE S ) 

+ (C 9 + Cio) {DE L + CEi)] , (9) 



AU 
^0 



= Ci 

+ c 2 



1 DEl + ^CP L + DA L + (DP L - DPl(c)) + (DPA L - DPA L {c)) 



2 -DE L - X -CE L + CA L + {CP L - CP L (c)) + (CPi L - CPA L {c)) 



(10) 



4 = Ci [DP L (c) + 5P4 L (c)] + C 2 [CP L (c) + CPi L (c)] 

+ C 3 [CP L + CA L + 2DA L + CP L + 2DP L + DP L (c) + DP L (s) 

+ CPA L + 2DPA L + DPA L (c) + DPA L (s)] 

+ C 4 [£>P L + 2CA L + DA L + 2CP L + CP L (c) + CP L (s) + DP L 

+ 2CP^ L + CP^l L (c) + CP^ L (s) + Z)P^ L ] 

+ C 5 [CE R + CA R + 2DA R + CP R + 2DP R + DP R (c) + DP R (s) 

+ CPA R + 2DPA R + DP4(c) + DPA R (s)] 

- 2C 6 [DEs + 2CAs + DA s + 2CPs + CPs(c) + CPs(s)+DPs (11) 
+ 2CPA S + CPA s (c) + CPi s (s) + DP A s ] 

+ hj 7 [CE R - DE R - CA R + DAr - CPr + DP R + 2DP R (c) - DP R {s) 

- CPAr + DPAr + 2£>PA fl (c) - DPAr(s)} 

+ C 8 [CP S - DE S - CA S + DA S - CP S - 2CP s (c) + CP s (s) + DP S 

- CPA S - 2CPA s (c) + CPi s (s) + DP A s ] 

1 By Fierz rearrangement we could have adopted a basis where only (V — A) x (V ± A) structures are 
present. We prefer however to work with operators written in terms of colour-singlet bilinears. 



+ ^C 9 [CE L - DE L - CA L + DA L - CP L + DP L + 2DP L {c) - DP L {s) 

- CPA L + DPA L + 2DPA L (c)-DPA L (s)] 

- ^Cio [CE L - DE L - CA L + DA L - CP L - 2CP L (c) + CP L {s) + DP L 

- CPA L -2CPA L {c) + CPA L {s) + DPA L ] . 



In eqs. (g)-(ll) we have not shown explicitly the argument \i of the Wilson coefficients 
Ci,...,Cio. Notice that also the diagrams are /i-dependent, since they correspond 
to contractions of renormalized operators Qi{ij). For the fields we have assumed the 
standard convention B + = bit, = bd, ir + = ud, tt~ = —du and tt° = —\j\[2(uu — 
dd). 

In Aq, penguin diagrams always appear in the GIM combinations DP^ — DPl(c), 
DPAl—DPAl(c), . . . (in the following we denote these combinations as GIM-penguins). 
Had we taken m t <C My/ and a renormalization scale [i larger than the top quark mass, 
i.e. mt <C /U <C Mw, similar GIM combinations would have appeared in Aq, namely 
DPi,(t) — DPl(c), etc. Since the physical value of mt is so large, the diagrams DPzit), 
CPi(t), etc. are replaced by complicated structures. These arise from the contractions 
of the penguin and electro-penguin operators Q3-Q10, originated in H e ff when we 
remove the top quark from the effective theory. In the literature, Aq only is identified 
as penguin pollution. We want to stress again that in the effective theory there are 
"penguin operators" Q3-Q10, which originate from the imperfect GIM cancellation oc- 
curring when n <S mt and "penguin diagrams" which arise from the Wick contractions 
of all the operators of T~t e ff- 

The coefficients of the penguin operators Q3~Qio are of order a s /127rln(m|/^ 2 ) 
and contain the short distance contribution (from scales between fi and mt) of the 
virtual top and charm quarks. For /i ~ m&, these coefficients are rather small, e.g. the 
dominant term is due to Qq, for which Cq/C2 ~ —0.03. Thus, unless the corresponding 
matrix elements are very large, the penguin operators are not expected to give large 
corrections, at least in the — ► tt + tt^ case. In the B® — ► 7r°7r° case, the relative 
correction due to penguin operators, Qq in particular, may be more important due to 
the large cancellations present in the amplitude (the so called "colour suppression" to 
be discussed below). 

This is not the end of the story, however. In A there are other terms, specifically 
those in the first line of eq. (|IT|). These terms, denoted as "charming penguins", come 
from penguin contractions of the operators Q\ and Q2, and have to be understood 
as long distance contributions in the matrix elements of these operators. Because of 
the unitarity relation A c = — X u — At, they give contributions to both Aq, in the GIM 
combination DPl-DPl(c), and to A . Unlike the case of the penguin operators Q3- 
Q10, the coefficients of Q\ and Q2 are of 0(1) and, since there is no reason to believe 
the corresponding matrix elements to be small 0, they are potentially relevant even for 
the B^ — > 7r + 7r~ amplitude. On the other hand We notice that their contribution can 
be further enhanced by the factor |r| ~ 2. 

Penguin diagrams are also present in Aq and can alter in A~* and A 00 the relative 
size of the term proportional to X u with respect to that proportional to A(. This can 



2 This is true in spite of the fact that, for (V — A) x (V — A) operators, penguin diagrams vanish if we 
assume factorization, as demonstrated by the important role that they are expected to play in kaon decays. 



be seen by looking to the expression of A introduced in eq. (Q) in terms of the different 
amplitudes 

A qA + ~ qX* u l + r*(A^-A^)/{A^-A§) 

P A+- p\ u i + T(Al-Ai)/(A%-A%) ' 1 j 



where 



q = VtdVtb = 1- p-iy = e - 2 
V V td V tb I- p + ir] 

K V ubVud = P~ i7 l = g -2i 7 

K v udV*b P + ir i 

r = J-P-^ (13 ) 
p + irj 



p and 7] being the CKM parameters in the Wolfenstein parametrization |16| and /3+7 = 
7r — a. 

At this point we have all the formulae needed for the study of the uncertainties on 
sin 2a and in the calculation of R. 

3 B { \ K°K°, B+ tt + K° and B Q d K n 

We now show that the B% K°K°, B + -> tt + K° and B° d -> K + ir~ amplitudes are 
dominated by GIM- and charming-penguin diagrams. Let us consider the three cases 
separately. For simplicity, all the formulae in this section are computed in the 577(3) 
symmetric limit, e.g. DPl{s) = DP^. 

B° d -> K°K° We start by studying this case for which all the relevant quantities 
have already been defined. The amplitude is dominated by GIM and charming-penguin 
diagrams because emission and annihilation diagrams are only produced by the inser- 
tion of the operators Q3-Q10 which have very small Wilson coefficients. By defining 
A{B° d -> K°K°) = X U G F /V2 x A K , the complete expression of Ax is given by 

A K = Ci [(DP L - DP L {c)) + tDP l (c) + (DP A - DPA L {c)) + rDPA L (c)] 

+ C 2 [(CP L - CP L {c)) + tCP l (c) + (CPA - CPA L (c)) + rCPA L (c)] 

+ r {C 3 [CE L + CA L + 2DA L + CP L + 3DP L + DP L (c) 

+ CPA L + 3DPA L + DPA L (c)] 

+ d [DE L + 2CA L + DA L + 3CP L + CP L (c) + DP L 

+ 3CPA L + CPA L (c) + DPA L ] 

+ C 5 [CE R + CA R + 2DA R + CP R + 3DP R + DP R (c) 

+ CPAr + 3DPAr + DPA r (c)\ 

- 2C 6 [DE S + 2CA S + DA S + 3CP S + CP s (c) + DP S 

+ 3CPA S + CPA s (c) + DP A s ] (14) 

+ ^C 7 1-CEr - CA R - 2DA R - CP R + 2DP R (c) 

- CPAr + 2DPA r (c)) 

+ C 8 [DE s + 2CA s + DA s -2CP s (c) + DP s 



- 2CPA s (c) + DP A s ) 

+ ^C 9 [-CE L -CA L -2DA L -CP L + 2DP L (c) 

- CPA L + 2DPA L (c)] 

- ^C w [DE L + 2CA L + DA L -2CP L (c)+DP L 

- 2CPA L (c)+DPA L ]} . 

By neglecting all penguin-operator contributions we get 

A K = C\ [(DP L - DP L (c)) + tDP l {c) + {DP A - DPA L (c)) + tDPA l (c)\ 

+ C 2 [(CP L - CP L {c)) + tCP l (c) + (CPA - CPA L {c)) + tCPA l (c)] .(15) 

This shows that unless the matrix elements of penguin operators, Qq for example, are 
exceedingly large (much larger than their estimates in the factorization hypothesis [3]), 
the amplitude is dominated by GIM- and charming-penguin diagrams. Equation (15|) 
is only given to display the relevant terms; in all our numerical calculations we always 
used the complete expressions for all the amplitudes. 

B+ _> Tr+if In order to study B + -» ir + K° and B° d -» K + 7r we have to in- 
troduce new quantities. For these decays, the operator Qf is replaced by Q[ u = 
(bs) fy-A\(uu)ty—A) with a Wilson coefficient C[ = C±, and similarly for all the other 
operators in (y): this comes from the fact that we are considering now AB = —AS = 1 
transitions instead of the AB = AD = 1 ones of sec. We also define A^ = V us V* b 
and similarly X' c and X' t : X' u is of 0(A ) whereas A^, and A^ are of 0(A 2 ), where A ~ 0.22 
is the sine of the Cabibbo angle in the Wolfenstein approximation. With the exception 
o£ A' uct , we will omit the superscript ' for the rest of this section. We write 

Gf 



A(B+ -> vr+K ) = -|(AU+ + + \' t A%) , (16) 



where 



A\ = C 1 [CA L + DP L ] + C 2 [DA L + CP L ] 
A\ = C 1 DP L (c) + C 2 CP L {c) 



A\ = -C 3 [CE L + CA L + CP L + 3DP L + DP L (c)} 

- C A [DE L + DA L + 3CP L + CP L {c)+DP L ] 

- C 5 [CE R + +CA R + CP R + 3DP R + DP R (c)] 

+ 2C 6 [DE s + DA s + 3CP s + CPs(c) + +DP s } (17) 

- ic 7 [-CE R + 2CA R - CP R + 2DP R (c)\ 

- C 8 [DE S - 2D A s - 2CP s (c) + DP S ] 

- ^C 9 [-CE L + 2CA L -CP L + 2DP L (c)} 

- ^Cio \-DE L + 2DA L + 2CP L (c) - DP L ] . 



By neglecting all penguin-operator contributions and Cabibbo suppressed terms we get 



A{B + -» tt + K») = -j=X c A° + = -^X' c [CtDPLic) + C 2 CP L (c)] , (18) 



By -> i^+vr" We write 



A(B° d - tf+7r> = ^ (AUo + AUo + , (19) 



where 



A* = C X [-CE L -DP L \ + C 2 [-DE L -CP L \ 
Al = -C x DP L {c)-C 2 CP L {c) 

4 = -C 3 [-CE L -CA L -CP L -3DP L -DP L (c)} 

- C i [-DE L -DA L -3CP L -CP L {c)-DP L ] 

- C 5 [-CE R - CA R - CP R - 3DP R - DPr(c)} 

+ 2C 6 [-DE s -DA s -3CPs-CP s (c)-DP s } (20) 

- ~C 7 [-2CEr + CA R + CP R - 2DP R (c)} 

- C s [2DE S - DA S + 2CP s (c) - DP S ] 

- X -C % [-2CE L + CA L + CP L - 2DP L (c)] 

- ^Cio [~2DE L + DA L - 2CP L (c) + DP L ] . 

By neglecting all penguin-operator contributions and Cabibbo suppressed terms we get 
A(B° d - K+7T-) = ^X' c Al = -^A' c [C^Px(c) + C 2 CP L (c)) , (21) 

In A(B + —>■ tt + K°) and A(B d — > i^ + 7r~), if charming-penguin diagrams are very small 
then both Cabibbo suppressed contributions and terms due to penguin operators must 
be included in the calculation since they are of the same size. 

4 Estimates of the diagrams 

In this section, we discuss the criteria adopted to evaluate the diagrams appearing 
in eqs. (|8|)-(|ll]) and in sec. ||. Notice that the value of the diagrams is given only 
once that the renormalization prescription (RP) and renormalization scale fi of the 
operators have been fixed. Under a change of RP or \i the values of the diagrams must 
be changed in such a way as to compensate the corresponding changes of the Wilson 
coefficients Cj, thus giving the same physical predictions up to and including next-to- 
leading logarithmic corrections 17]— 1|22[| Q The state of the art in the calculation of the 



3 Notice that under a change of the RP or of /i, contributions attributed to matrix elements of some 
operators can go into the Wilson coefficients or the matrix elements of others, and viceversa. 



matrix elements of the operators is such that, given the complexity of the expressions 
in eqs. (jg|)-(|ll|), this turns out to be impossible. For example, factorized amplitudes 
are RP and scale independent, being expressed in terms of physical quantities. Thus 
they cannot compensate a variation of the coefficient functions. On the other hand, 



the possibility that lattice calculations be able to compute (g)— (11) with sufficient 
accuracy appears to be rather remote. This is particularly true for A 00 where delicate 
cancellations are likely to occur between different contributions, see also the discussion 
and the end of this section and section |5[ In the following, in order to check the 
stability of the results, at fixed values of the diagrams we vary the Wilson coefficients 
by changing RP and by taking 2 GeV < /i < 10 GeV. 

We now discuss the assumptions made in the evaluations of the different diagrams: 

1) DP A and CPA These are Zweig suppressed diagrams which we assume to give 
a negligible contribution. 

2) Electro-penguins In order to monitor the effects of the electro-penguins, we only 
consider the contributions coming from the operators Qg and Qio since these 
operators have coefficients much larger than Qj and Qs- 

3) DEl and CEl In most of the theoretical analyses, these diagrams give the largest 
contribution to the B^ — > tttt amplitudes. If only emissions are present, there 
are three independent quantities namely \DEl\, \CEl\ and arg(DEi x CE* L ). 
Without loss of generality we can then write CE L = £,DE L e i5 t. We vary 0.0 < 
£ < 0.5: this range covers the value preferred by the analysis of D-meson two- 
body non-leptonic decays, which suggests £ ~ 0, and the value derived from 0,2/ a\ 
extracted from B —* Dit and B — > Dp decays [||, 23]. Moreover it includes the 
canonical value £ = 1/N C where N c is the number of colours. On the basis of 
some estimates of FSI |24|, [25|] , we do not expect the relative phase 5^ to be larger 
than ~ 0.5, at the energy scales flowing into the pion system in U-decays. We 
enlarged the range of 5^ up to ~ 1 in order to check which kind of effect could be 
produced by a large phase. 

4) DAl and CAl These are diagrams usually neglected since arguments can be 
made that, in the factorization hypothesis, they are suppressed by a factor Jb/Mb 



(besides colour suppression in CAl) [15]. Jb/Mb is related to the 5-meson wave 



function in the origin. A further suppression factor comes from the matrix element 
of the divergence of the vector current which creates the pion pair. In ref. pfifl , 
however, a fit to two-body .D-meson decays, resulted in a non negligible value 
for the annihilation diagrams, corresponding to DAl/DEl ~ 0.3(m<j — m^)/ f n . 
It is not clear how these results should be scaled to the £>-meson case and for 
degenerate quark masses. 

Rescattering effects, which have been shown to persist even for large quark masses 
[24], can also enhance the value of annihilation diagrams with respect to factor- 



ization estimates [27], see also [28|. As discussed in ref. [27], emission and anni- 
hilation diagrams are connected via FSI. For example, the CA-diagram can be 
seen as a DE followed by rescattering of the final states. In the 1/N C expansion, 

since the scattering amplitude is of order 1/N C , this gives immediately the cor- 

1/2 

rect leading dependence on the number of colours for CA^. In fact DEl ~ N c , 

— 1/2 

CAl ~ N c and CAl ~ t/aDEl, where t\a is proportional to the scattering 



amplitude, which is of 0(1/N C ). The only potentially large contribution from an- 
nihilation diagrams comes indeed from the term oc C2CAL in eq. (^), since all 
the other annihilation contributions (either DAl or CAl) have (much) smaller 
coefficients. Moreover, several arguments can be made to show that the value of 
DAl is expected to be at most of the size of CAl and its largest contribution 
multiplies C\ which is about 1/5 of C2. In our numerical analysis, we have ex- 
plicitly checked, by varying DAl between zero and CAl-, that its effect is rather 
marginal. For this reason, in sec. |5| we only discuss the case with DAl = 0. To 
take into account rescattering effects, we parametrize CAl as CAl = iDEl^a- 
t]a is a complex "inelasticity" coefficient, the absolute value of which has been es- 



timated to be of the order of some tenth [ 24 , ^] . In our numerical study we take it 
real with < t\a < 0.5 The same rescattering mechanism relates disconnected 
and connected diagrams in the penguin case, which is discussed below. 

5) GIM-penguins Penguin- like contractions appear both in ^4q and A . Here follow- 
ing we discuss the two cases separately. In Aq we always find the combinations 
DPl — DPl(c) and CPl — CPl{c) and one could argue that, because of the 
large final state energy at disposal, GIM cancellation makes these contributions 
negligible. The GIM mechanism, however, is expected to be effective only at 
short distances, i.e. when a high momentum flows in the penguin loop. For low 
momenta, i.e. if we look to long-distance effects, these diagrams can also be 
interpreted as emission diagrams followed by rescattering. For example, 

CP L - CP L (c) ~ DE L {B° d -» ^7r)5 (vr^ -> tttt) 

-DE L {B Q d -^DD)S (DD -> vrvr) 
~ Uq»(n(pB-q)\JZ b \B° d )S (™^<Kn) 

-f D q»(D(p B - q)\J^ b \B° d )S (DD - tttt) , (22) 

where So is the strong interaction S'-matrix and Jjf , ' cb are the weak vector cur- 
rents. Since fv/ f-n ~ 1-5 [§9| and the form factors relative to (tt(pb — Q)\JjT\B d ) 
are expected to be smaller than those relative to (D(pb — Q)\J^ c \B d ) ||-||, 
by about a factor of 2, F cb (M|,) ~ 2F ub (M%), it is not clear how effective is 
the GIM cancellation between the two contributions. A large cancellation be- 
tween the charm and up contributions may still take place if the relative factor 
Sq(D D — ► 7T7r)/5o(7r7r — > tttt) compensates the differences due to phase space, 
decay constants and form factors [J 

The discussion of rescattering effects for GIM-penguins strictly follows that of 
the annihilation diagrams made in 4). For this reason, given our ignorance of So, 
we use the parametrization CPl — CPl(c) = iDElt/p, with < r]p < 0.5, and 
ignore the contribution from DPl — DPl(c). 

6) Other penguin diagrams In A , penguin contributions are not GIM suppressed. 
Thus we expect that rescattering effects play a minor role. Penguin contributions 
are of two kinds: either they correspond to the insertion of left-left operators, Q%, 



4 



All estimates give Re t]a 3> Im r\A W% p5[ 



5 In this respect a combined measurement of the B — > tttt and B — > D D branching ratios would be very 
interesting. 



Q2, Q3, Qi, Q9 and Q10, or they are given by the insertion of Q5 and Qq, for 
example in DPr or DP$. 

In the previous section, we have noticed that penguin contractions of the operators 
Qi and Q2 can give large effects since the corresponding coefficients are of O(l). 
In our numerical analysis we find that a modest relative phase between DPl and 
CPl can have dramatic effects. For this reason we introduce two parameters r/x, 
and 5 L , by writing DP L = \DE l \vl and CP L = e iSL \CE L \vL, with 0.0 < f] L < 1.2 
and < 5l < 0.5. The range of values of t]l is dictated from the fact that there 
is no reason to expect very large/small values for these matrix elements, so we 
take them of the same order than the corresponding emission diagrams. As for 
the phase, we limit the maximum of 5l to 0.5, as we did for 5%. For simplicity, we 
take r]L flavour independent, i.e. we use the same value of t]l for DPl, DPl(s) 
and DPl(c). We do not connect CPl to DEl, as we did for GIM-penguins, 
because, as said above, we do not have to advocate long-distance effects coming 
from emissions followed by rescattering. In general, we should consider a complex 
value of tjl- We checked, however, that the largest effects come from the relative 
phase between DPl and CPl and for this reason, only the case with tjl real will 
be discussed in the following. 

Left-right penguin contributions only appear together with all the other contrac- 
tions of Q5 and Qq. The latter have the same topology of the diagrams considered 
so far, but different chiral structures. Rather than introducing another set of free 
parameters for the (V — A) x (V + A) and (S + P) x (S — P) diagrams, we consider 
globally the matrix elements of the operators Q5 and Q§ and write 

(7nv\Q 5 \B° d ) = rj 5 CE L , {^\QM) = ^DE L . (23) 

In kaon decays, there is a common prejudice that the operators Q5 and Qq trig- 
ger the octet enhancement 6 . If really the explanation of the enhancement relies 
on the matrix elements of these operators (as also suggested by lattice calcula- 
tions |30[j ), then 775,6 can be as large as 5. Since the kinematical configuration 
is so different in S-decays, due to the large mass of the 6-quark, and many of 
the arguments for the enhancement are based on the chiral expansion, we do not 
really expect 775,6 to be as large as in the kaon case. Since we cannot exclude, 
however, values of O(l), we vary 0.0 < 775,6 < 2.0. 

We summarize the discussion above, by giving the expressions (|8|)-(|i"T|) in units of 
DEl, and in terms of the parameters £, <5^, 6l and of the r/^s 

A\ = _I(C 1 + C 2 )(l + ^) (24) 
A\ = -I(C 9 + C 10 )(l + £e W *) , (25) 

3 3 s 

6 In the kaon case, if \x is larger than the charm mass, which plays the role of the top quark mass for GIM 
effects, the operators Q§ and Qe are hidden in the matrix elements of Q\ — Q\ and Q2 — Q2 (with b — > s in 
the operators of the effective Hamiltonian) . 
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Before discussing the numerical results of our analysis, we want to add some ob- 
servations about colour-suppressed processes. The starting point are eqs. (j8|)-(|i~l|) 
[fl24|)— (p7|)1, which contain the expressions of the relevant amplitudes. For the sake of 
illustration, let us consider first the case where only DEl and CEl are non vanishing 
and we neglect penguin operator contributions. In this case, when we insert (|S|)-(|iT|) 
in eqs. (g) and @ the following combinations occur 



oc CxCEl + C 2 DE L = C 2 DE L ( 1 + 
A w oc C 1 DE L + C 2 CE L = C 1 DE L (l + ^e iS i 



(28) 
(29) 



Numerically, the ratio —C\/C 2 ~ 0.2-0.3 is approximatively equal in size and opposite 
in sign to the expected value of £. This implies that the second term in eq. fl28| ) is a 
small correction (of the order of 10%) to the first one while the two terms in eq. ( |29| ) 
tend to cancel (for small values of the phase <5g). This, together with the smallness of 
Ci, is at the origin of colour suppression. In D-decays, where a very small value of £ 
is preferred, colour suppression is expected to be less effective f\. 

As can be read from eqs. d^)-(p|), a similar colour-suppression pattern is present 
in all the AI = 1/2 annihilation and penguin diagrams generated by the operators 
Q\-Q 2 , which are the operators with the largest coefficients: 



00 



oc d[DA L + (DP L 
+ C 2 [CA L + {CP L - 



DP l (c))+tDP l (c)] 
CP l (c))+tCP l (c)} + 



(30) 



Thus, if all the connected diagrams are about a factor of £ smaller than the correspond- 
ing disconnected ones, it is very difficult to resurrect the — > 7r°7r° amplitude, unless 
the matrix elements of the penguin operators Q3-Q10 are very large. The same argu- 
ments, however, could be applied to kaon decays, with the surprising conclusion that 
\Aq I ~ 2 1 A 2 1 in sharp disagreement with the experimental observation \Aq\ S> \A 2 \. It 
is then clear that dynamical effects play a fundamental role in this game. 

Of course, it is possible to argue that the dynamics of the decay for kaons and 
i?-mesons is completely different and that factorization is a very bad approximation 
for light mesons. It is clear, however, that one cannot exclude some residual dynamical 



7 This could alternatively mean that factorization is a poor approximation. 



effect which enhances the B® — ► 7r°7r° amplitude over its factorized value. The contri- 
bution of Qi and Qi to Aq, non-zero phases 5^ and di or the connected-disconnected 
diagram-exchange mechanism for CAl and CPl — CPl{c) discussed above may provide 
such effects. 



5 Numerical results 

In this section, we present numerical results for the quantities sin 2a and R, that only 
depend on ratios of amplitudes. The latter are all proportional, via the parameters £, 
. . ., 776, to DEl, the value of which cancels out in the ratios. 

One could, however, also be interested to know the variation of BR(B® — > ir + ir~) for 
different assumptions about GIM-penguin, charming-penguin, annihilation diagrams 

etc. Thus we also give below = BR(B® — > it + it~) / BR(B® — > 7t + 7t~)|.d.e l , where 

BR(B® — > 7r + 7r~) is the branching ratio computed for a given set of the parameters 
£, . . ., 77q, whilst BR(B® — ► 7r + 7r~)|£>£ L is the branching ratio with all the diagrams 
put to zero, but DE^. In this way, the reader can use his preferred model to compute 
D£x and predict the physical value of BR(B® — > 7r + 7r~). In the following, R will also 
denote BR(B% -> tt°tt )/BR(B° -> 7r+7r-)| DBi . 

For the determination of a and the relative error, we define 

sin2ae // = ^, (31) 
A 

which is the quantity that can be extracted from the time-dependent asymmetry by 
measuring the coefficients of cos(AM^t) and of sin(AMrft) in eq. ([l]). The uncertainty 
on the "true" value of sin 2a can be estimated by constructing 

A = sin 2a e ff - sin 2a . (32) 

Finally we present our results for the ratios R\ = BR(B® — > K°K°)/BR(B® — > 
^ + 7r-)|D£ L , ^2 = -► ir + K°)/BR{B° -► tt+tt-)!^ and # 3 = Bfl(B|j -► 

K + tt~)/ BR(B® — > 7r + 7T _ )|£)E £ , the values of which have a strong dependence on the 
contribution of charming-penguin diagrams. 



In subsec. |5.1| the results for R and i2_| under different assumptions on the operator 



matrix elements are presented; in subsec. |T3 the effects of these assumptions on the 



determination of sin 2a are analyzed, while subsec. |5.3| contains the numerical results 
for Ri, i?2 and R3. 



5.1 Estimates of R and R^ 

Given the large number of parameters, in order to understand which diagrams give 
important/unimportant contributions, for several values of £ and 5% we vary one single 
of the ?7jS (and 5l) of the set defined in sec. |I| at the time, while keeping all the others 
fixed to zero. We have also checked that the main features of the results discussed 
below are the same irrespectively of the choice of the RP and of the renormalization 
scale n (2 GeV < /x < 10 GeV). 
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Ci/C 2 


p (GeV) 


LO 


NDR 


HV 


RI-Landau Gauge 


2.0 


-0.31 


-0.27 


-0.32 


-0.24 


5.0 


-0.19 


-0.16 


-0.19 


-0.15 


10.0 


-0.12 


-0.09 


-0.12 


-0.08 






p (GeV) 


LO 


NDR 


HV 


RI-Landau Gauge 


2.0 


0.52(7) 


0.49(8) 


0.53(7) 


0.46(8) 


5.0 


0.42(8) 


0.39(9) 


0.42(9) 


0.39(9) 


10.0 


0.36(9) 


0.33(9) 


0.36(9) 


0.32(9) 



Table 1: Values ojC\jCi ond£ for different RP and at different values of the renormalization 
scale p. The values of £ have been computed assuming 02/01 = 0.25 ± 0.05, from B — > Dtt, 
Dp, D*tt and D*p decays, and using £ = (02/01 — C\/C>i)/{\ — 02/01 x C1/C2). 



a) Colour suppression for DEl and CEl Colour suppression depends on the value 
of the Wilson coefficients, £ and the phase <5g. Information on the value of £ is 
usually obtained by analyzing several -B-decay channels (B — * Dtt, Dp, D*tt and 
D*p |23[| ) and is correlated to the Wilson coefficients used in the analysis. In 
most of the cases, the leading order (LO) coefficients at p = mj are used and 
a value of £ ~ 0.4 is found. From a comparison of the LO results, with those 
obtained at the next-to-leading order (NLO) in different renormalization schemes 
(NDR, HV and RI p2]]), we find that the value of £ is not very sensitive to the 
(considered) RP, whereas it can vary from 0.23, at p = 10 GeV, to 0.60, at p = 2 
GeV 0. The correlation between the ratio Ci/C 2 and the extracted value of £ is 
shown in table |l|. We have taken 02/01 = 0.25 ± 0.05, in agreement with most 
of the analyses of the experimental data [23]. In absence of a fully consistent 



treatment of the amplitudes at the NLO, which could be obtained if some lattice 
calculation existed p3| , we have no reason to prefer a particular value of the scale, 
or to maintain any correlation between p and the value of £. For this reason, in 
the following we take the coefficients computed at the LO for p = 5 GeV from 
ref. j22| but vary £ in the range 0.2 < £ < 0.5, which covers almost all the values 
m table |. 

We have a comment on the choice of the Wilson coefficients made in some recent 
studies which may be useful to the reader. In refs. |5|, [H]] and [^] they used 
the renormalization scheme independent coefficients introduced in ref. pPR , and 
computed for the full basis (Q) in ref. [31]. Though perfectly legitimate, this 



choice corresponds to a value of the ratio Ci/C 2 = —0.27 (and to £ ~ 0.47), 
sensibly larger than those found at leading order or at the NLO in the NDR, HV 
or RI schemes at a renormalization scale p ~ m&, see for example ref. [22] Q. 



On the basis of the discussion of the previous section, see eq. (|29|) , for 5^ = 



8 This determination of £ is only indicative, since it is not known whether the same value of £ should be 
used for B — > Dir and B — * tttt decays. 

9 Indeed the renormalization scheme of ref. jl9| has never been completely specified, because the external 
states on which the renormalization conditions were imposed were not given explicitly in the paper. 
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Figure 3: The ratio R is given as a function of £ at three values of 5%. In the amplitudes, 
only contributions of emission diagrams are included. 

the ratio R has a minimum for £ ~ —C1/C2 ~ 0.2 and rises for larger values of 
£. The suppression is softened for non- vanishing values of 5t, as can be seen in 
fig. |3|, and a value of R = 0.03-0.06 is obtained for £ ~ 0.4-0.5 and 5^ = 0.5 (to 

be compared to R = 0.02 at £ = 0). Notice that only varies by less than 20 

% in the range of £ and 5^ considered here (it would only vary by about 10 %, as 
a function of 5%, at £ = 0.4 fixed). 

b) Annihilation and GIM-penguin diagrams The effect of annihilation diagrams de- 
pends on 5^ and £. Let us first fix £ = 0.4 and vary 5^. For small values of 6%, up to 
5^ ~ 0.4, R increases with x\a and reaches a value of ~ 0.1 at t]a = 0.5; for larger 
values of 5%, instead, there is destructive interference which balances the positive 
effect of 5^ on colour suppression. The situation is illustrated in fig. |j. Contrary 

to R which is subject to appreciable variations, with the same parameters R-\ 

changes by less 20 % even at the largest value of r\A- Similar behaviours, but on 
different ranges of 5^ are observed at other values of £. 

The effect of GIM-penguin diagrams is essentially the same as for annihilation 
diagrams, since, as can be immediately seen from eq. (p6|), the largest Wilson 
coefficient C2 multiplies the combination tja + np. 
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Figure 4: The ratio R is given as a function of tja at three values of 5^ . In the amplitudes, 
only contributions of annihilation and emission diagrams are included. 

c) Q§ and Qq and Electro-penguins Q$ and Qe give relatively small corrections both 

to R and , because their coefficients are very small. Of the two terms, the 

contribution of Q5, which has the smaller coefficient and is colour suppressed, 
is always very small. Q§ has the effect of changing by about ~ 80% and ~ 5% 

the ratios R and respectively when tjq = 2. A large effect from Qq can 

only be obtained at extremely large values of its matrix element. In comparison, 
electro-penguin diagrams always give tiny corrections. 

d) Charming penguins The most important of all the effects is given by charming 
penguins. The explanation was already given in the previous section: these dia- 
grams correspond to the insertion of the operators Q\ and Q2 which have large 
Wilson coefficients and the contribution of which is also augmented by the fac- 
tor r. They may easily enhance the — ► amplitude, and change also 
appreciably the — > 7r + 7r~ rate. 

In figs. |H and [| we give R and -R_| as a function of ni at three value of 6l 

(for £ = 0.4 and 5^ = 0). We observe that, for t]l ~ 1-0, there is a quite 
substantial effect when 5l is different from zero: R can reach values as large 
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Figure 5: The ratio R is given as a function ofr/i at for 5l = 0, 0.25 and 0.50, at ^ = 0.4 and 
5^ = 0. In the amplitudes, only contributions of emission diagrams and charming penguins 
are included. 



as 0.25 (for 5 L = 0.25) and even 0.35 for 5 L = 0.5 Even with S L = we 
find a large enhancement, which can lead to R = 0.10-0.15 at relatively modest 

values of t\l- In this case there is even an increase of about 20% of This 

discussion shows that values of BR(B^ — ► 7r 7r°) as large as 1-3 x 10 -6 (assuming 
BR(Bj — ► 7r + 7r~) = 1-2 x 10 -5 ) are indeed easily possible. 

5.2 Determination of sin 2a 

In this subsection, we discuss the uncertainties in the determination of sin 2a. These 
uncertainties are parametrized in terms of the shift A introduced before. A is computed 
at different values of sin 2a which are obtained by varying the CKM weak phase 5. The 
values of of sin 2a are computed from the expression sin 2a = Im(r/T*), with r given 
in the last of eqs. (|l3|), and using p = a cos 5 and n = a sin 5 with a = 0.36 and 
—0.8 < cos 5 < +0.7 (sin 5 > 0). In the figures below, we give A as a function of 

10 The effect is striking for the B® — > ir + ir~ rate which at i]l — 1-2 and Sl ~ 1.0 almost vanishes. We 
believe, however, that this is a quite remote possibility. 
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Figure 6: The ratio _R + _ is given as a function ofrji, forSi = 0, 0.25 and 0.50, at£ = 0.4 and 
8% = 0. In the amplitudes, only contributions of emission diagrams and charming penguins 
are included. 



sin2a e " because the latter is the quantity which is measured experimentally. The 
strange behaviour of A for sin2a e -^ k5 0.80 comes from the fact that two different 
values of cos 5, for cos 5 ^ — 0.6, correspond to the same value of sin 2a. Had we 
varied cos 5 within the range allowed by the combined analysis of the K°-K° and of 
the B®-B® mixing amplitudes [ 33 1 , i.e. —0.3 < cos 5 < 0.9, the two-fold ambiguity 
would have disappeared, because the interval in cos<5 limits, in this case, sin 2a to 
values smaller than about 0.9. 

In the following, we discuss together the ) of subsec. lO] and the case d) 

separately. 

a)-c) In fig. [7| we show a band of possible values for A. This band has been obtained 
in the following way: a) by varying £ between 0.0 and 0.5 with 5^ = 0, 0.5 and 
1.0, and all the other rji's taken to be zero; b) by fixing £ = 0.4 and varying 
rjA between 0.0 and 0.5 with 5^ = 0, 0.5 and 1.0 (np = 775 = rje = rji = 0); c) 
by taking £ = 0.4, <5 ? = and 0.0 < r/ 6 < 2.0 (va = np = 775 = ni = 0). We 
see that for all the cases a)-c) the size of A is about —(0.05-0.15) in almost the 
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Figure 7: A band of possible values for A as a function of sin 2a is given. This band 
corresponds to different choices of the parameters introduced in l)-5), see sec. [|. 



whole interval of values of sin 2a This happens because the terms depending 
on £, 5%, r\A and 775,6 give small variations (of the order of 20% at most) to the 
— * ir + ir~ amplitude. These results are similar to those found in refs. H 32] 
where different hypotheses or approximations where used. 



d) We varied 771, and 5l as at the point d) of subsection |5,1| , while taking £ = 0.4 
and and all the other 77jS to be zero. In this case the uncertainty on sin 2a 
can be very large, corresponding to A ~ 0.4-0.8, see fig. |8[ The reason, as 
already discussed above, is that charming penguins give large corrections to the 
— ► 7r + 7r - amplitude too, while in the other cases considered above [a)-c)], the 
corrections are large only for the B^ —> amplitude where large cancellations 
occur. 

In order to guide the reader, we end this subsection by giving in table [2] a set of 
numerical results for sin 2a e ^ and A obtained for different choices of the parameters 
introduced in a)— d) (also reported in the table). In the set of examples given in table 
A £ 0.6. 



Smaller values are only found when sin la is close to +1 or — 1. 
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Figure 8: A family of curves for A is given as a function of sin 2a. They correspond to 
different choices of the parameters r\ L (r\ L = 0, 0.5, 1.2) and 5l (Sl = 0, 0.25, 0.5) introduced 
in d), see sec. f|. 



5.3 Charming-penguin dominated branching ratios 

Our main results for the ratios of rates introduced before are summarized in table 
||, where a large set of possibilities has been considered. We have introduced a label 
"No." to distinguish the different cases which will be discussed below 0. The table is 
instructive because it shows the strong correlation between R and -Ri,2,3> particularly 
when charming-penguin contributions are large. We will make use of this correlation, 
combined to the experimental information BR(B® — > K + ir~) = (1.5^q'| ± 0.2) x 
10~ 5 and BR(B® — > n + ir~) < 1.5 x 10 -5 H, and to prejudices based on theoretical 
estimates of the matrix elements, in order to reduce the uncertainties. Given the 
present precision in the measurement of BR(B® — > K + ir~), and our ignorance on 
the matrix elements, the discussion below is, at the moment, preliminary and semi- 
quantitative and only intended to illustrate a method which may become very useful 
when the experimental precision will increase and theoretical predictions for GIM- and 
charming-penguin diagrams will become available. 

12 In order to check the stability of the results, other cases have been analyzed but not given in table 
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0.24 


8 


0.32 


0.35 


-0.56 


0.21 


9 


0.24 


0.27 


-0.64 


0.13 


14 


0.17 


0.21 


-0.67 


0.11 


15 


0.53 


0.57 


-0.36 


0.42 


18 


0.21 


0.25 


-0.63 


0.14 


23 


0.52 


0.56 


-0.38 


0.39 


27 


0.43 


0.47 


-0.47 


0.3 



Table 2: We give results for sin 2a e ^ and A obtained with different choices of the parameters 
5^, . . ., rjiS. In all cases £ = 0.4. In the table, the value of cos 5 used in the different examples, 
and the corresponding value of sin a, are also shown. The label "No. " , which is the same of 
table 0, allows a study of the correlation between A and the values of the branching ratios 
discussed in subsec. [57^ . 

We proceed as follows: 

i) BR(B® — > 7r + 7r~) Different theoretical estimates of this branching ratio, obtained 
by using the factorization hypothesis but with different models to evaluate the ma- 
trix elements, are consistent within a factor of two and give BR(B® — > 7r + 7r~) = 
1-2 x 10~ 5 ||-||]. The spread of these predictions is mostly due to differences in 
the values of the form factors of the weak-currents and in the value of £ used for the 
calculation of the emission diagrams. The results, instead, are only marginally 
affected by the contribution of the penguin operators Q3-Q10, which in some 
cases have been omitted (charming-penguin diagrams have never been consid- 
ered). This happens for the following reason: in the factorized case the matrix 
elements of Q5 and Qq, which are expected to give the largest contributions of 
all penguin operators, correspond to 775^ ~ 1 for which the correction is of a 
few percent, as can been seen by comparing No. 2 to No. 5 in table ^. This 
explains the relative stability of the theoretical results. We then assume that 
1-2 x 10 -5 is a range representative of the theoretical uncertainty for emission 
diagrams 0. On the basis of the above discussion, in the following we consider 

13 Predictions on BR(B { j — > 7r°7r°) are, instead, very unstable because in this case even a small error on 



three possibilities BR(B% -» Tr + 7r-)\ DE+C E = (0.5, 1.0, 1.5) x 1(T 5 . The lat- 
ter two cases have been used because they are in the ballpark of all estimates, 
the first value because it accounts for large deviations from the factorization 
predictions. Values of BR(B d — > tt + tt~)\de+ce larger than 2 x 10~ 5 are ex- 
cluded on the basis of ii) below. BR(B d — ► tt + tt~)\de+ce corresponds to the 

case No. 2 of table ^, where = 0.85. Thus, for example, if we assume 

BR(B° d -» 7r + 7r-)|ujg+c£ = 1.5 x 10" 5 , the ratio R = 0.54 (No. 18) gives 
BR{B° d -» 7r°vr ) = 0.54/0.85 x 1.5 x 10~ 5 = 0.95 x 10~ 5 . 

ii) BR(B° d -* K+-K-) vs BR(B G d -> tt+tt') From CLEO we expect BR(B% -► K+ir-) 
to be larger than BR(B d — ► 7r + 7r _ ). Thus we exclude all those cases where 
Si?(5^ -» K+7T-) < BR(B d -> vr+TT-), namely No. 1-6, 13, 17, 18, 21, 22, 25, 
26. 

iii) BR{B° d -» K+7T - ) We also exclude all cases where 5^(5^ -> ET+vr^) < 1 x 10" 5 
or BR(B° d -> K+tt~) > 2 x 10~ 5 . 

Indeed, from the experimental results, BR(B d — > K + ir~) ~ 2.5 x 10~ 5 or £#(#2 -> 
i^+vr") ~ 0.5 x 10~ 5 and also BR(B° d -> iT^) ~ -> tt+tt') are still open 

possibilities. However, we want to stretch here the experimental constraints in order 
to show how this analysis works. 

We now discuss the results of table ||, where the main results of the selection 
based on i)— iii) are given. For comparison, in this table we also give the results for 
No. 2, where only emission diagrams contribute. We first notice that the charming- 
penguin enhancement of BR(B + -> ir + K°) and BR(B% -> K+ ir ) is always very 
large, corresponding to a factor of 10-40, with respect to standard case No. 2. A 
stable prediction of the results in the table is that we also expect BR(B + — > ir + K°) 
to be comparable to BR(B d — > /T + 7r - ) and of the order 1 x 10~ 5 . Although, on the 
basis of table ^, a much large enhancement of BR(B d — > 7r°7r°) would be possible, 
the constraints i)— iii) limit the effect of charming penguins so that the final result is 
about a factor of 2 larger than previous estimates ||, 0] . A larger enhancement is still 
possible, however, if we relax the selection constraints. For example, without ii), case 
No. 18 is acceptable and, for BR(B d — ► n + ir~)\DE+CE = 0-5 X 10~ 5 , it gives 



BR(B° d 


— ► 7T + 7T 


) = 1.2 Xl0~ 5 , 


5#(5° vr°vr ) = 0.3 x 10~ 5 , 


BR(B° d 




) = 0.3 x 10~ 5 , 


Si?(5+ -> vr+i^ ) = 0.6 x 10~ 5 , 


BR(B° d - 




) = 0.8 x 10~ 5 . 


(33) 


These results g 


;ive a lar£ 


«e BR{B° d -> 7r°7r°), 


and are still compatible with the experi- 



mental measurement of BR(B d — > X + 7r ) and the present limit on BR(B d — ► 7r + 7r ). 



This example also show the importance of the reduction of the experimental error for 
determining charming-penguin effects. 

We finally notice that for BR(B d — ► X°^°), besides charming penguins, GIM- 
penguins are also relatively important. Moreover GIM- and charming-penguin dia- 
grams appear in the r.h.s. of eq. ( [To]) in the same combination as in Aq, eq. (p7|). A 

the evaluation of a single diagram is amplified by the large cancellations of the different contributions. This 
is also shown by comparing different entries in the table, where variations of R over one order of magnitude 
are observed. 



measurement of BR(B® — ► K°K°) is then very important for the determination of the 
combined value of GIM and charming-penguin diagrams in B® — > tttt decays and may 
lead to a reduction of the uncertainties in the determination of sin 2a. Unfortunately, 
in most of the cases, this branching ratio remains rather small, BR(B° d -► K°K°) ~ 
5 x 10~ 7 , in spite of the enhancement due to charming penguins. For this reason this 
may be a difficult measurement even for S-factories. However, from the example dis- 
cussed in flU and the results of the table, a BR(B% -> K°K°) = 2-3 x 1CT 6 is still 
compatible with the present experimental constraints and the experimental search of 
this decay mode is very important. 



6 Conclusions 

In this paper we have discussed several effects which could enhance the B% — ► 7r°7r° 
branching ratio. Among these, we have shown the presence of diagrams involving 
operators containing charmed quarks, denoted as "charming penguins", which were 
never studied before. Since there is no reason a priori to expect the value of these 
diagrams to be small, and the corresponding Wilson coefficients are of 0(1), and not 
of 0(a s /127rln(m^/// 2 )), they may have a large effect in the B d — > 7r°7r° decay rate 
and in the determination of the weak phase a extracted from the measurement of the 
B® — > 7r + 7r~ time-dependent asymmetry. In absence of any realistic estimate of the 
value of charming-penguin diagrams, we allow them to vary within reasonable ranges 
of values and show that, without further constraints, they can lead to a BR(B d — > 
7r°7T ) ~ 1-3 x 10 6 . Correspondingly, the uncertainty in determination of sin 2a from 
the time-dependent asymmetry could be as large as 0.4-0.8. 

We have also shown that GIM- and charming-penguin diagrams dominate the 
B d — > K°K° amplitude because emission and annihilation diagrams are only produced 
by the insertion of the operators Qs~Q io which have rather small Wilson coefficients. 
GIM- and charming-penguin diagrams appear in the B° d -> K°K° amplitude in the 
same combination as in the B d — > tttt case. Thus, if we assume SU(S) symmetry, a 
measurement of BR(B d — > K°K°) contains important information on these contri- 
butions and may help to reduce the uncertainty in the extraction of sin 2a from the 
B® — ► 7r + 7r _ time-dependent asymmetry. 

A striking effect of charming-penguins occurs in B + — > ir + K° and B d — > K + tt~ 
decays. In this case, the enhancement easily gives values of BR(B + — ► tt + K°) and 
BR(B d — ► K + tt~) of about 1 x 10~ 5 , 10-40 times larger than those obtained by 
considering emission diagrams only. This is because in the latter case the main con- 
tribution is Cabibbo suppressed. Our findings are supported by the recent results of 
the CLEO Collaboration Q. Theoretical predictions of the charming-penguin (and 
GIM-penguin) amplitudes, either from the lattice, or with QCD sum rules, or with any 
other non-perturbative technique, are then urgently needed. 
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Table 3: We give results for R = BR(B° d -> ir ir°)/BR(B% -> 7r+7r-)| DBi; 7? + _ = 
- n + n-)/BR{B° d - 7r+7r-)| DEi; i? x = - K°K )/BR(B° d - tt+tt-)!^, 

i? 2 = -> n + K°)/BR{B° d -> 7T+7r-)| DSi and i? 3 = -> K+tt-)/ BR(B° d -> 

7r+7r_ )|DE L ; obtained with different choices of the parameters . . ., rjiS. 



No. 


5^ -» vr+vr- 




Bl K°i^ 


B + _> vr+X u 


S^ -► K+TT- 




Bi? x 10 5 


Si? x 10 5 


Bi? x 10 5 


Si? x 10 5 


BR x 10 5 
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Table 4: Branching ratios for B° d -> tt+tt-, £g -> 7r°7r° ; 5^ -> i^if , fi+ -> vr+iT and 
5° — > K + 7i~ . The cases are labeled as in table [|. The BRs are normalized assuming three 
different values of BR(B d — > 7r + 7r~) \de+ce- 



